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SUMMARY 

Attitude and heading determination using GPS interferometry is a well-understood 
concept. However, efforts have been concentrated mainly in the development of robust 
algorithms and applications for low dynamic, rigid platforms (e.g. shipboard). This paper 
presents results of what is believed by the authors to be the first realtime flight test of a 
GPS attitude and heading determination system. The system is installed in Ohio 

University’s Douglas DC-3 research aircraft. Signals from four antennas are processed by 
an Ashtech 3DF 24-channel GPS receiver. Data from the receiver are sent to a 
microcomputer for storage and further computations. Attitude and heading data are sent 
to a second computer for display on a software generated artificial horizon. Demonstration 
of this technique proves its candidacy for augmentation of aircraft state estimation for flight 
control and navigation as well as for numerous other applications. 


INTRODUCTION 

The Global Positioning System (GPS) will soon be capable of providing worldwide 
positioning and timing services to properly equipped users. In addition, interferometric 
techniques can be used to obtain realtime, highly accurate attitude and heading information 
from the GPS measurement data as well. GPS interferometry is a technique based on the 
measurement of phase differences of the GPS carrier signal observed at two or more 
antennas. Attitude and heading of a platform can be calculated from measurements with 
two linearly independent interferometers where each interferometer consists of two separate 
antennas. 

GPS interferometry has been used extensively in precise static and kinematic 
positioning. Experimental results of static relative positioning using GPS interferometry 
were first reported in 1982 (ref. 1), showing subcentimeter accuracies for baseline lengths 
between 8 and 52 meters. Further stationary experiments were presented in 1983 (ref. 2). 
After that time, implementation issues were emphasized, such as the development of fast 
initialization algorithms (ref. 3). Realtime low dynamic shipboard tests were performed in 
1988 (ref. 4), with reported angular accuracies on the order of several degrees for a 60- 
centimeter baseline. The first non-realtime determination of aircraft pitch and heading 
using one GPS interferometer was reported in 1989 (ref. 5), revealing pitch and heading 
uncertainties below 1 milliradian for a 23-meter baseline. Today, sub-milliradian accuracies, 
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even for 1 -meter baselines, are attainable for rigid platforms using off-the-shelf equipment 
(ref. 6). In this paper, we present methods and results for realtime pitch, roll, and heading 
determination of a dynamic, flexible platform. 


PRINCIPLES OF ATTITUDE DETERMINATION 

In the formalism of rigid body kinematics, attitude is defined by the rotation transformation 
which relates a coordinate system fixed in space to a coordinate system fixed in the body 
(ref. 7). For the purposes of aircraft navigation, the space coordinate system typically is 
defined to be locally level. The positive x-axis is pointed towards north, the positive y-axis 
is pointed east and the positive z-axis is pointed down. Finally, since the transformation is 
concerned with rotation only (and not translation), the origin of the space coordinate system 
is coincident with the aircraft center of gravity. Similarly, the body (or aircraft) coordinate 
system also has its origin at the aircraft center of gravity. However, for the aircraft 
coordinate system, the positive x-axis is directed along the fuselage axis towards the nose, 
the positive y-axis points along the right wing (right according to the pilot), and the positive 
z-axis is directed towards the floor of the aircraft (ref. 8). 

The two coordinate systems are related through the antenna locations. Under the rigid body 
assumption, the antenna locations remain fixed in the aircraft coordinate system. Their 
locations in this system are determined prior to the flight through standard GPS static 
surveying techniques. While in flight, the antenna locations in the space coordinate system 
are determined using GPS interferometry. Thus at any given epoch, the antenna locations 
in the two coordinate systems are related as follows: 

[R\ s =mm a (i) 

where: [R] s is the matrix containing the locations of the antennas in space coordinates. [R] a 
is the matrix containing the locations of the antennas in aircraft coordinates. [T] is the 
rotation transformation matrix. 

The transformation is a 3x3 matrix containing the nine direction cosine elements. It is the 
product of three rotation matrices each of which specifies a rotation about a certain aircraft 
coordinate system axis. Thus, 

m=[<H[0]['jf] ( 2 ) 

where: [<j>] represents a rotation about the fuselage axis (roll); [0] represents a rotation 
about the wing axis (pitch); and [ijr] represents a rotation about the vertical axis (heading). 

The transformation matrix may now be expressed in terms of the three rotations (also 
known as the Euler angles): 
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Determination of the instantaneous values of the Euler angles proceeds as follows. First, 
the matrix containing the antenna locations in aircraft coordinates is formed prior to the 
start of attitude determination. Secondly, realtime GPS interferometric measurements are 
used to form the matrix containing the antenna locations in the space coordinate system. 
Solving for the transformation matrix involves postmultiplying the space coordinate matrix 
by the generalized inverse of the aircraft coordinate matrix: 




(4) 


Finally, the Euler angles can be calculated as follows: 

e=sin- 1 (-r 31 ); 4>=sin- 1 (^); t=sin 1 (-%) (5) 

COS0 COS0 

The calculations resulting in this Euler angle triad have been presented to provide insight 
into the attitude determination process. However, although this triad is probably the best 
understood means of attitude and heading expression, it is far from robust and arguably 
should not be considered for use in a realtime integrated navigation system (ref. 9). This 
lack of robustness can be seen most easily by examining equations (3) and (5) for the case 
of 0 = ±90°. Note that although 0 can be recovered, 0 and f are lost. As a result, the 
Euler angle triad is not an all-attitude parameterization. Although many parameter! zati on s 
of [T] are possible, quaternions offer conciseness, efficiency, and stability not possible with 
other parameterizations (refs. 8-9). 

Further practical considerations involve the aircraft structure. The derivation given above 
rests on the assumption that the antenna locations are fixed in the aircraft coordinate 
system. Error mechanisms are thereby introduced by the non-rigid body behavior of the 
aircraft. The fuselage by itself is a good rigid body approximation, but inclusion of the wings 
(with their normal flexing) make the overall aircraft far from being so. An additional 
consideration is more subtle. The derivation assumes that the location of the aircraft center 
of gravity is known relative to the antenna positions. Although a nominal location can be 
used, the true location will change over the course of the flight. 


GPS INTERFEROMETRY 

Consider a simplified GPS interferometer, as depicted in figure 1, consisting of two GPS 
antennas and the necessary electronics to determine the instantaneous phase of a received 
GPS carrier signal. Usually, the distance between the antennas and a GPS satellite is much 
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larger than the baseline length between the two antennas, such that the GPS carrier signal 
can be treated as a plane wave. Under the assumption that the satellite clock and the 
receiver clock are synchronized, the phase difference A 0 between simultaneous carrier phase 
measurements at the two antennas can be written as 

A 0 = 0j - 0 2 = b *2. + NX (6) 

where 0 j and <p 2 are the measured carrier phases in wavelengths for antennas 1 and 2 , 
respectively; b is the three-dimensional baseline vector pointing from antenna 1 to antenna 
2; 6 is the unit vector pointing from the interferometer to the satellite; and N is an integer 
number accounting for the ambiguity due to the limited length of the carrier wavelength k. 
The phase difference between antennas is also referred to as a Single Difference (SD) (ref. 
10). The determination of the unit vector to the satellite requires the approximate 
coordinates of the interferometer and the satellite coordinates. The latter are calculated 
from the Navigation Data transmitted by the satellites. A standard GPS navigation solution 
is used to determine the interferometer coordinates. The baseline solution is not very 
sensitive to the interferometer coordinates, as was shown by Jurgens (ref. 11), because of 
the large distance to the satellite with respect to the relatively short baseline length of the 
interferometer. 

If the satellite clock and the receiver clock are not synchronized, each of the phase 
measurements, (p l and 0 2 , contain an unknown satellite clock phase offset and a receiver 
clock phase offset. If the phase measurements are made simultaneously, then the phase 
difference A 0 does not contain the unknown satellite clock phase offset. Actually, since the 
time of transmission for each of the phases transmitted by the satellite will be slightly 
different by an amount At=(A 0 A.)/c, a small offset will be present in the phase difference 
due to the change in the satellite clock offset during the interval At. This offset can be 
corrected using the satellite clock coefficients transmitted by the satellite in the Navigation 
Data. 

The unknown receiver clock offset can be eliminated by taking the difference between two 
independent Single Differences, also referred to as a Double Difference (DD) 

DD = A0j - A0 2 = ^•(i i -l 2 ) + (A,-A 2 )A. (7) 

Thus a Double Difference is a linear combination of four phase measurements obtained 
from two satellites at two antennas, as shown in figure 2. This can be extended to eliminate 
the unknown integer ambiguity number by differencing two Double Differences at two 
epochs to obtain a Triple Difference. The Triple Difference is useful for static positioning 
but only provides the rate of change of the Double Difference for dynamic applications. As 
such it does not provide an accurate instantaneous solution for the baseline vector. 

Assuming that the integer ambiguities are known, then using equation ( 6 ), it follows that 
three independent Double Differences are required to solve for the baseline vector b = 

(b 1 ,b 2 ,b 3 ) T : 
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Several techniques have been proposed to quickly determine the integer ambiguity for each 
of the Double Differences (refs. 12-13). In general, search methods are very effective to 
quickly resolve the integer ambiguities whenever redundant measurements are available. 
The basic idea is to evaluate the measurement residuals for potential baseline vector 
solutions. The techniques used are very similar to those applied to GPS fault detection (see 
ref. 14 for a summary of these techniques). If, for instance, signals from five satellites are 
received by one GPS interferometer, four independent Double Differences can be formed. 
Equation (8) is rewritten as: 

DD = Hb (10) 


where DD is a 4-by-l Double Difference measurement vector and H is a 4-by-3 data matrix 
containing the satellite unit vector differences. The residuals are the differences between 
the actual Double Difference measurements and those predicted based on the least squares 
solution for b, the approximate interferometer coordinates, and the known satellite 
coordinates. The least squares residual vector, e r , is obtained by projecting DD onto a space 
perpendicular to the estimation space spanned by the columns of H. Thus, e r is orthogonal 
to the estimation space. An efficient method to obtain e r is by decomposing the H matrix 
into the product of a real orthonormal matrix Q and an upper triangular matrix R using the 
QR factorization (refs. 15-16). Replacing H in equation (10) by the product QR and pre- 
multiplying both sides by Q T yields (Q T Q = I) 

R b = Q t DD (11) 

The rank of R is equal to the rank of H; therefore, the bottom row of R consists of zeros 
only. Equation (11) can be divided into two equations: 
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( 12 ) 


U b = Q l DD or b = U~ l O x DD 


0 = a DP 


where U consists of the first three rows of R, Qj consists of the first three rows of Q T , and 
g is the transpose of the fourth row of Q T . The top expression of equation (12) relates the 
measurements to the baseline coordinates, forming the least squares solution. The vector 
<j and the Double Difference vector DD are othogonal. The latter property is used to 
obtain the norm of the residual vector if measurement errors are present on the Double 
Differences: 


/ \ 
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The vector q has to be calculated only once for a search at one particular epoch. Thus, to 
calculate the residual for a potential solution, only four multiplications and three additions 
are required. 

Some prior knowledge is required about the interferometer in order to reduce the potential 
number of solutions for p. Several methods can be used, including: 

• Limit the baseline length of the interferometer to, for instance, one meter such that 
the integer ambiguity number must lie in the interval [-5,5] (one GPS wavelength, k, 
is approximately 0.19 meters); 

• Use a priori information about the orientation of the interferometer such that the 
search can be limited to correspond to the approximate baseline orientation. 
Approximate attitude and heading information could for instance be obtained from 
a low cost three-axis accelerometer and magnetometer; 

• Use Double Differences calculated from the GPS pseudorange measurements (code 
phase Double Differences) to obtain approximate values for the carrier phase Double 
Differences. 

If the Double Difference uncertainty can be limited to ±5 A, then even an exhaustive search 
would only require the calculation of (ll) 4 = 14,641 residuals, which would require seven 
times that many or 102,487 floating point calculations. This can easily be accomplished 
within one second using a microcomputer. In this paper, we estimate the initial DD 
ambiguities by filtering the relatively noisy code phase Double Differences with the carrier 
phase Double Differences using a complementary Kalman filter (see also ref. 17). In the 
complementary filter, the carrier phase DDs are used to provide a reference for the code 
phase DDs. The performance of the complementary filter is very similar to the Hatch filter 
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(ref. 18), but has much lower computational and storage requirements. Although only one 
DD is estimated, the complementary filter model can easily be extended to include the rate 
of change of the DD as well. The equations for the (scalar) complementary Kalman filter 
are summarized below 


£ = £-1 + (PD k - DD k _j) 
Pk = Pk-\ + Q 
K = Pk (Pk + r )~' 


& = 4 + K( \ ~ £ ) 


Pk = (1 - *k)A 


where & denotes the estimated code phase DD with variance q; r is the variance of the code 
phase DD measurement z; and p is the variance of the estimation error. The first two lines 
in equation (14) propagate the estimated code phase DD and the estimation error variance 
from time (k-1) to time k, followed by the calculation of the Kalman gain (line 3), and the 
update of the estimated code phase DD and the estimation error variance (lines 4 and 5). 
Figure 3 shows an example of both code phase DDs and carrier phase DDs for a 20-minute 
section of a flight test. The DDs were obtained from a GPS interferometer on the fuselage 
of a DC-3 research aircraft during roll maneuvers. The next section provides a detailed 
description of the aircraft test bed. Note that the code phase DDs and the carrier phase 
DDs are separated by an unknown integer wavelength ambiguity. Figure 4 shows both the 
carrier phase DDs and the filtered code phase DDs. As expected, the filtered code phase 
DDs follow the dynamics very well. Finally, the difference between the carrier phase DDs 
and the filtered code phase DDs are shown in figure 5. After the filter is initialized with 
the first code phase DD, it takes typically less than 20 seconds for the carrier phase DD to 
be within ±5X with respect to the filtered code phase DD. The remaining differences 
between the filtered code phase DDs and the carrier phase DDs are due to the code phase 
multipath from the aircraft structure. 

The filtered code phase DD is used to calculate an initial estimate for the absolute carrier 
phase DD. Next, an exhaustive search is performed to find those combinations of integer 
ambiguities that will result in low residuals. Depending on the carrier phase measurement 
noise and the satellite measurement geometry, more than one potential set of ambiguities 
may be found. Further eliminations of potential sets are achieved by applying known 
baseline length constraints. For an interferometer mounted on a semi-rigid platform, the 
baseline length can be determined to within a few millimeters. The knowledge of the norm 
of the baseline vector is limited by differential antenna phase center migrations, differential 
multipath and structural flexing of the platform. Next, additional redundant measurement 
data can be used to resolve the integer ambiguities, again by examining the residuals. 
Although statistically unlikely, it would still be possible to select the wrong set of 
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ambiguities; however, by continuing to examine all measurement residuals throughout the 
experiment, the resulting inconsistencies in the data would be quickly detected. 

Finally, after the integer ambiguities have been determined for each of the Double 
Differences, they can be used to maintain aircraft attitude and heading for long periods of 
time until loss-of-lock occurs. It is also important to add that once the integer ambiguities 
are resolved, redundant measurements are no longer necessary to maintain attitude and 
heading. Assuming known baseline lengths, a rigid platform would only require continuous 
lock on the signals from three satellites at three antennas. In this case, to find the 
approximate interferometer coordinates, another sensor input, such as the baro-alti meter 
or a stable frequency reference source, is required. If the baselines are subjected to 
structural flexing, continuous lock on four satellites at three antennas is required. 


DESCRIPTION OF THE TESTBED 

The realtime attitude and heading determination testbed consists of three main components: 
1) Aircraft, antennas, and cabling; 2) GPS receiver; 3) Data storage, computation, and 
display units, as depicted in figure 6. Four microstrip antennas have been installed on Ohio 
University’s Douglas DC-3 research aircraft, N7AP. One antenna is located on each wingtip 
and tv/o antennas are located on top of the fuselage with one near the nose and the other 
towards the tail (fig. 7). RG-400 cable is used because of its high phase stability and small 
physical diameter. Phase stability is essential since each of the four cables is approximately 
25 meters in length. The cables feed the signal from the antenna preamplifiers to the GPS 
receiver located in a standard equipment rack installed in the cabin midway between the 
nose and tail of the aircraft. 

The GPS receiver being used is the Ashtech Three-dimensional Direction Finding (3DF) 
24-channel receiver (ref. 6). The 3DF is configured as four 6-channel sections with each 
section receiving signals from a separate antenna. The receiver is packaged in a ruggedized 
3-5/8" X 8-1/2" X 8" housing and requires a 12 watt power source. Realtime measurement 
data at a selectable rate of up to once per second (both raw and processed) are available 
over two serial ports at data rates up to 115,200 baud. 

The storage, computation and display units consist of two microcomputers. A 25 MHz 
80386-based computer receives realtime data from the 3DF, performs storage, computation 
and sends realtime attitude information to a second computer. The second computer is a 
12 MFIz 80286-based PC. It generates an artificial horizon and heading indicator and drives 
it with the attitude information supplied by the 80386 machine. 

Baselines between the antennas were determined using standard GPS static surveying 
techniques. The lengths of these baselines are given in Table 1 . These measurements have 
also been verified using a steel tape measure and were found to agree to within 5 

millimeters. 
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Baseline 

Length in meters 

ant 1 - ant 2 

7.781 

ant 1 - ant 3 

14.468 

ant 1 - ant 4 

14.466 

ant 2 - ant 3 

14.941 

ant 2 - ant 4 

14.924 

ant 3 - ant 4 

28.194 


Table 1. Antenna baseline lengths for 
Douglas DC-3 test bed. 


GROUND and FLIGHT TESTING 

Ground and flight testing of the GPS attitude and heading determination system started in 
April of 1991. To date, three demonstration flights have been conducted (April 18, May 8 
and May 22). Each flight test lasted for approximately 1.5 to 2 hours, and was preceded by 
a data collection session on the ramp at Ohio University’s airport (UNI) in Albany, Ohio. 
During the first two flights only raw measurement data were collected which were used for 
the continued development of realtime attitude and heading determination algorithms. 
During the third flight, long periods of realtime attitude and heading determination were 
achieved. Moreover, the 3DF GPS receiver maintained continuous lock on at least 3 
satellites at 3 antennas during all maneuvers. This is an important result, because it allows 
for continuous attitude and heading determination. We also found that whenever the 
receiver would loose lock on a particular satellite due to temporary shielding by the 
airframe, the receiver would typically re-acquire the signal within a few seconds. 

Aircraft flexing was found to be the main reason for discontinuities in the attitude and 
heading data. The carrier phase DD residuals would increase which caused the algorithm 
to re-initiate the integer ambiguity search. This could be avoided by allowing larger DD 
residuals in combination with the use of a higher baseline length uncertainty. 

A representative section of the ramp test results collected from the 3DF receiver on May 
22, 1991, is shown in figures 8 through 10. It should be noted that the aircraft was not static 
during this data collection session; small changes in pitch, roll and heading angles were 
introduced by wind gusts and by the flight crew moving around inside the aircraft. In 
addition, some multipath effects can be observed in the data. During this session, five 
healthy GPS satellites were visible (satellites 2, 6, 13, 14 and 15). The Position Dilution of 
Precision (PDOP) was approximately equal to 3 throughout the test. The noise level on the 
roll and heading angles is on the order of 0.1 mrad (0.006 degrees). Pitch angle noise is 
approximately three time larger because of the shorter baseline used for the pitch angle 
determination. Absolute accuracy results are not available due to the lack of a reference 
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for the attitude data which must be aligned with the phase centers of the GPS antennas. 
However, previous comparisons between the 3DF and kinematic GPS results were found 
to agree to within 0.1 mrad (ref. 6). 

Realtime data collected from the 3DF receiver during two sections of the third flight test 
conducted in the evening of May 22, 1991 are presented next. The noise levels on the pitch, 
roll and heading angles are difficult to analyze directly, because of the many attitude and 
heading variations with relatively large magnitudes due to normal flight. Instead, the 
emphasis of the results presented in this section is on the capability of GPS interferometry 
to provide realtime, in flight attitude and heading information. Figure 11 shows the ground 
track of a 27-minute section of the third flight which was conducted in the vicinity of Ohio 
University. Figures 12 through 14 show the aircraft heading, roll and pitch angles as a 
function of time. Note that the noise-like behaviour of the angles are all due to actual 
aircraft dynamics (the period of the angle variations is on the order of 10 to 20 seconds). 
Approximately 12 minutes after the start of the run, attitude and heading information was 
lost for approximately 30 seconds. As discussed earlier, this was caused by poor DD 
residuals which caused the attitude and heading algorithm to re-initiate the integer 
ambiguity search. Shown in figure 15 is the lock history of the 3DF receiver; even in the 
presence of signal shielding due to airplane maneuvers, the receiver maintained lock on a 
minimum of four satellites at 3 antennas. 

Figures 16 through 20 show the results for a 13-minute section of the third flight test during 
which a series of 10 and 30-degree roll maneuvers were performed. Continuous attitude and 
heading information was available using four or five GPS satellites. During the 30-degree 
rolls, only one of the five visible satellites was not observed continually by at least three 
antennas. 


SUMMARY, DEVELOPMENTS AND APPLICATIONS 

The feasibility of using realtime GPS interferometric attitude and heading determination for 
dynamic semi-rigid platforms has been demonstrated. During the course of this research, 
several issues have been identified which require further study. Among these are continued 
development of robust fault detection and isolation algorithms for the detection and 
correction of cycle slips, characterization and mitigation of multipath, and integration with 
an inertial measurement unit. 

GPS interferometry has numerous applications including: reliable aircraft/ spacecraft state 
estimation; measurement of structural flexing of large space structures; synthetic aperture 
radar (ref. 5); aerial photogram metry; pointing systems; and rapid alignment of inertial 
navigation systems for submarines and aircraft while maneuvering. 
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Figure 3. Code phase (bottom) and carrier phase (top) Double Differences. 



Figure 4. Carrier phase Double Differences (top) and filtered code phase Double 
Differences (bottom). 
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Figure 5. Difference between carrier phase Double Differences and filtered code phase 
Double Differences. 
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Figure 6. Block diagram of realtime attitude and heading determination test bed. 
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Figure 7. Douglas DC-3 GPS antenna locations. 



Figure 8. Ramp test heading residuals. 
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Figure 9. Ramp test roll residuals. 



Figure 10. Ramp test pitch residuals. 




Figure 11. Ground track for a 27-minute section of the third flight test. 


FLIGHT TEST: HEADING 



Figure 12. Flight test heading as a function of time. 




Figure 14. Flight test pitch as a function of time. 
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Figure 15. Flight test lock history. 



Figure 16. Ground track for a 13-minute section of the third flight test. 




FLIGHT TEST: HEADING 



Figure 17. Flight test heading as a function of time. 



Figure 18. Flight test roll as a function of time. 




Figure 19. Flight test pitch as a function of time. 
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Figure 20. Flight test lock history. 




